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We analyze the properties of a class of improved lattice topological charge density operators, constructed by 
a smearing-like procedure. By optimizing the choice of the parameters introduced in their definition, we find 
operators having (i) a much better statistical behavior as estimators of the topological charge density on the 
lattice, i.e. much less noisy; (ii) a multiplicative renormalization much closer to one; (iii) a large suppression of 
the perturbative tail in the corresponding lattice topological susceptibility. 



1. Introduction. 

In QCD an important role is played by the to- 
pological properties. By the axial anomaly, ma- 
trix elements or correlation functions involving 
the topological charge density operator q{x) can 
be related to relevant quantities of the hadronic 
phenomenology. We mention the topological su- 
sceptibility X) which is determinant in the solu- 
tion of the U(1)a problem, and the on-shell nu- 
cleon matrix element ofq(x), which can be related 
to the so-called spin content of the nucleon. 

Lattice techniques represent our best source of 
non-perturbative calculations, however investiga- 
ting the topological properties of QCD on the lat- 
tice is a non-trivial task. In a lattice theory the 
field is defined on a discretized set and therefore 
the associated topological properties are strictly 
trivial. One relies on the fact that the physical 
continuum topological properties should be reco- 
vered in the continuum limit. 

Considering a lattice version of q{x), q^{x), 
the classical continuum limit must be in gene- 
ral corrected by including a renormalization func- 
tion. In pure QCD, where q{x) is renormalization 
group invariant. 
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where Zi^g^) is a finite function of the bare cou- 
pling (jfQj going to one in the limit ^ Q. But at 
3q ~ 1, where simulations are usually performed, 
it may be very different from one. The finite re- 
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(n^,y(a;) is the product of link variables ?7^(a;) aro- 
und a 1 X 1 plaquette) is indeed very small at 
gl ~ 1: for SU{i) Z{gl = 1) ~ 0.18. 

The relation of the zero-momentum correlation 
of two q^ operators, Xl > with the topological su- 
sceptibility X is further complicated by an un- 
physical background term, which eventually be- 
comes dominant in the continuum limit. Indeed 
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Neglecting terms 0(a®), the background term 
Mi^g^) can be written in terms of mixings with the 
unity operator (so-called perturbative tail sca- 
ling as ~ a°) and with the trace of the energy- 
momentum (scaling as ~ a'*). In the case of 
the operator (2) and for 5?7(3), Mi^g^) is alre- 
ady dominant at (/q ~ 1 (it is about 85% of at 
3o = !)• As a consequence, using the heating me- 
thod to evaluate Zi^g^) and Mi^g^) [1], the uncer- 
tainty on X can be hardly pushed beyond ~ 10%. 

Another problem, which has come up in some 
studies concerning the lattice determination of 
the on-shell proton matrix element of q{x) [2], 
is that the lattice operator (2) is very noisy: a 
very large statistic and therefore expensive simu- 
lations are necessary in order to get a reasonable 
uncertainty on the final result. In view of a full 
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QCD lattice calculation the search for a better 
estimator appears a necessary step. 

We study, within the field theoretical approach, 
the possibility of improving the lattice estimator 
of q{x) with respect to all the problems listed 
above, that is we look for local lattice versions 
of q(x) which are less noisy, have a multiplicative 
renormalization closer to one, and whose corre- 
sponding Xl iio^ dominated by the unphysical 
backgroung signal M((jfo) in the region (jfg — 1- 
(Any Xl defined from a local will eventually 
be dominated by its perturbative tail in the conti- 
nuum limit. To the purpose of evaluate x a good 
result would be to have it small at (jfg — Ij which 
should be already in the scaling region.) 

2. Improved topological charge density 
operators. 

Inspired by the widely used smearing techni- 
ques, we consider the following set of operators 
defined in terms of smeared links V^^\x) 
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where n|i*2 is the product of smeared links V^^\x) 
around a 1 x 1 plaquette. Such smeared links are 
constructed by the following procedure: 
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where v[2{x) = vS'\x - 1^)1 V^'\x) and the- 
refore q^^\x) depends on the parameter c, which 
can be tuned to optimize the properties of q^^\ 
All these operators have the correct classical con- 
tinuum limit, i.e. for a ^ q^^\x) a'^q[x). 

Notice that the size of qi^\x) increases with 
increasing the integer parameter i, but they can 
still be considered as local operators when kee- 
ping i fixed while approaching the continuum li- 



mit. Anyway, as we shall see that a good impro- 
vement with respect to q^^\x) = q^ix) is already 
achieved for small values of i, by optimizing the 
choice of the parameter c. 

The procedure (5) may be used to improve any 
local operator involving link variables. Smearing 
methods to improve lattice estimators have been 
already widely employed in the study of long di- 
stance correlations, such as large Wilson loops 
and hadron source operators. 

3. Perturbative analysis. 

We have calculated Zi^g^) to one loop for the 
once-smeared operator q^^\x) within the Wilson 
action formulation. We find Z^^^g^) = 1 + 2130 + 
0{gt), where 
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At c = we recover the non-smeared result: Z = 
1 - 0.5353^ + 0{gl) ioi N = 2, Z = 1 - 0.908^^ + 
0(3q) for TV = 3 [3]. As c varies, the following 
extreme values of Z are obtained: zi = —0.136 
at c = 0.650 for TV = 2; zi = -0.247 at c = 0.677 
for TV = 3. In both cases, the reduction of zi is 
quite large, making the one loop estimate more 
reliable for typical values of g^ . 

For q^^\x), we have also calculated the lowest 
perturbative contribution to the mixing with the 
unity operator P^g^), which is a part of the ba- 
ckground term M^g^). 
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g 3TV(TV^ - 
^° 1287r4 
0.002867- 0.017685C+ 0.048665c^ 
0.075362c^ + 0.068526c^ 
0.034433c^ + 0.007445c^ 
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The minimum of this everywhere-concave polyno- 
mial is p{c = 0.872) = 1.4 X 10"^ Thus, for all 
TV, the leading order of P diminishes by more 
than two orders of magnitude compared to its 
non-smeared value (c = 0). 
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In the presence of dynamical fermions one 
should take into account the fact that, unlike 
pure gauge theory, the topological charge density 
mixes under renormalization with . The non- 
renormalizability property of the anomaly in the 
MS scheme means that the anomaly equation sho- 
uld take exactly the same form in terms of bare 
and renormalized quantities. However the renor- 
malization of dfij^[x) and q{x) is nontrivial [4]. 
A renormalization group analysis leads to the fol- 
lowing relation valid for all matrix elements of a 
lattice version q^{x) of q(x) in the chiral limit [5]: 

{i2Nfq,) = Y{gl){R) , (8) 

where Y^g^) is a finite function of g^, and 

(R) ^ {d,jlixU-)e.p r -^d-g (9) 

is a renormalization group invariant quantity, 
where dfij^[x)R— indicates the operator dfij^[x) 
renormalized in the MS scheme, and the func- 
tion j{g) is related to the anomalous dimension 
of the continuum operators q{x), df^j^^x) in the 

MS scheme: 7(3) = ^^Nfg^ + O [g^]. No- 
tice that (i?) is what can be naturally extracted 
also from experimental data. 

In perturbation theory one finds Y^g^) = 1 + 
(•^1 + J/i)fi'o + 0{gQ), where zi is the coefficient of 
the 0(30) term of the finite renormalization ofq^ 
in the pure gauge theory (cfr. Eq. (6)), and j/i 
turns out to be a small number: j/i = —0.0486 
for TV = 3 and Nf = 4. 

4. Non-perturbative analysis by the hea- 
ting method. 

Estimates of the multiplicative renormaliza- 
tions of the operators q^^\x) and of the ba- 
ckground term in the corresponding Xl t)e 
obtained by using the numerical heating me- 
thod [6,7], without any recourse to perturbation 
theory. We applied the heating method to the 
operators q^^\x) for i = 1,2 and for a number of 
values of c in the region < c < 1. We restric- 
ted our analysis to the SU{2) pure gauge theory, 
expecting no substantial differences for TV = 3. 
The measurements were performed at /3 = 2.6 



(3^ = 1.5384...), which is a typical value for the 
SU{2) simulations with the Wilson action. The 
estimates of Z^^\l3 = 2.6) from the plateaus ob- 
served when heating an instanton configuration 
are reported in Table 1, and should be compa- 
red with the value Z{l3 = 2.6) = 0.25(2) for the 
standard operator (2) [8]. 



Table 1 

2'(^)(/3 = 2.6) for i = 1,2 and some values of c. 



i 


c = 0.6 


c = 0.8 


c = 1.0 


1 


0.52(3) 


0.57(2) 


0.55(2) 


2 


0.67(2) 


0.74(2) 


0.68(2) 



This analysis confirms the one-loop perturba- 
tive calculations, that is the improved operators 
we considered have a multiplicative renormaliza- 
tion closer to one than that of the initial opera- 
tor q^{x). From Z{l3 = 2.6) ~ 0.25 of q^{x), we 
pass, optimizing with respect to the parameter c, 
to Z^^\(3 = 2.6) ~ 0.57 by one improving step, 
and Z^'^\(3 = 2.6) ~ 0.74 by two improving steps. 
For larger i we expect to get Z^'^^ closer and closer 
to one. On the other hand, we should not forget 
that increasing the number of improving steps the 
size of the operator q^^\x) increases, therefore one 
should find a reasonable compromise taking into 
account the size of the lattice one can afford in 
the simulations. 

A comparison of the above results for i = 1 
with the one-loop calculation shows that the con- 
tribution of the higher perturbative orders is still 
non-negligible, but not so relevant as in the case 
of the operator without improving. 

Another important property of the improved 
operators, we can infer from the heating method 
results, is that they are much less noisy than 
q^{x). A quantititive idea of this fact may come 
from the quantity e(^) = A2'(^)/2'(^), where AZ^'^ 
is the typical error of the data in the plateau du- 
ring the heating procedure. We indeed found for 
c ~ 1.0 and for an equal number of measurements: 
e(°)/e(i)~7and e^/ef^) ~ 18 . 
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An estimate of the background signal M^g^) 
can be obtained by measuring Xl '^^ ensembles 
of configuration constructed by heating the flat 
configuration [7,8]. The estimates of M^'\l3 = 
2.6) from the plateaus observed in the heating 
procedure are reported in Table 2, and should be 
compared with the value M(/3 = 2.6) = 2.09(5) x 
10~^ relative to the standard operator (2). 



Table 2 

M(^)(/3 = 2.6), multiplied by a factor 10^ for 
i = 1,2 and some values of c. 



i 


c = 0.6 


c = 0.8 


c = 1.0 


1 


0.59(2) 


0.36(1) 


0.26(1) 


2 


0.23(1) 


0.13(1) 


0.07(1) 



Notice the strong suppression of the backgro- 
und term in the zero-momentum two-point corre- 
lations constructed with the improved operators. 
For c ~ 1 the reduction is about a factor 8 when 
performing one improving step, and about 30 by 
two improving steps. For larger values of i, the 
suppression is expected to be larger and larger. 

The suppression of the background term in 
Eq. (3) together with the relevant increase of Z 
should drastically change the relative weights of 
the contributions to Xl do — i-^- 
relevant region for Monte Carlo simulations. For 
example if the improvement for SU(3) is simi- 
lar to that achieved for SU{2), using the optimal 
operator for i = 2 at gg = 1 the unphysical term 
in Eq. (3) may become a small part of the total, 
allowing a precise determination of x by the field 
theoretical method. 

5. Conclusions. 

We have analyzed the properties of a class of 
improved lattice topological charge density opera- 
tors. Such improved operators look promising for 
the lattice calculation of the on-shell proton ma- 
trix element of the topological charge density ope- 
rator in full QCD, which is related to the so-called 



proton spin content. Indeed their use should over- 
come the difficulty due to the large noise obser- 
ved in preliminary quenched studies [2] , and their 
multiplicative renormalizations are much closer to 
one. 

They should also provide a relevant impro- 
vement in the determination of the topological 
susceptibility by the field theoretical method in 
the SU(3) gauge theory, since the unphysical ba- 
ckground term should get strongly reduces while 
the term containing x gets enhanced by the larger 
values of the multiplicative renormalization. This 
should allow a precise and independent check 
of the alternative cooling method determinations 
(see e.g. Ref. [9]), whose systematic errors are 
not completely controlled. Furthermore the im- 
proved operators may open also the road to a 
more reliable lattice investigation of the behavior 
of the topological susceptibility at the deconfine- 
ment transition, where cooling technique does not 
seem to give satisfactory results [10]. 
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